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Abstract 

In cellular automata with memory, the unchanged maps of the conven- 
tional cellular automata are applied to cells endowed with memory of their 
past states in some specified interval. We implement Rule 30 automata 
with a majority memory and show that using the memory function we 
can transform quasi-chaotic dynamics of classical Rule 30 into domains of 
travelling structures with predictable behaviour. We analyse morphologi- 
cal complexity of the automata and classify dynamics of gliders (particle, 
self-localizations) in memory-enriched Rule 30. We provide formal ways of 
encoding and classifying glider dynamics using de Bruijn diagrams, soliton 
reactions and quasi-chemical representations. 

Keywords: Rule 30, memory, elementary cellular automata, chaos, com- 
plexity dynamics, gliders, emergency. 

1 Introduction 

Elementary cellular automaton (CA) is a one-dimensional array of finite au- 
tomata, each automaton takes two states and updates its state in discrete time 
depending on its own state and states of its two closest neighbours, all cell up- 
dates their state synchronously. A general classification of elementary CA was 
introduced in [41]: 

Class I. CA evolving to a homogeneous state. 



Class II. CA evolving periodically. 
Class III. CA evolving chaotically. 

Class IV. Include all previous cases, also know as class of com- 
plex rules. 

In this classification Class IV is of particular interest because rules of the 
class exhibit non-trivial behaviour with rich and diverse patterns, e.g. Rule 54 
CA [23, 24]. 

2 Basic notation 

2.1 One-dimensional cellular automata 

One-dimensional CA are represented by an infinite array of cells xi where i £ Z 
and each x takes a value from a finite alphabet S. Thus, a sequence of cells {xi} 
of finite length n represents a string or global configuration conS. This way, the 
set of finite configurations will be represented as An evolution is represented 
by a sequence of configurations {cj} given by the mapping $ : S™ — > E n ; thus 
their global relation is provided as follows 

$(c 4 ) -> c t+1 (1) 

where t is time and every global state of c is defined by a sequence of cell states. 
Also the cell states in configuration c* are updated at the next configuration 
c* +1 simultaneously by a local function ip as follow 

(2) 

Wolfram represents one-dimensional CA with two parameters (k,r). Where 
k = |E| is the number of states, and r is radius of ncighboourhood. Elemen- 
tary CA are defined by parameters (k — 2,r = 1). There are S™ different 
neighborhoods (where n = 2r + 1) and k k different evolution rules. 

In computer experiments we were using automata with periodic boundary 
conditions. 

2.2 Cellular automata with memory 

Conventional cellular automata are ahistoric (memoryless): i.e., the new state 
of a cell depends on the neighborhood configuration solely at the preceding time 
step of ip (eq. 2). 

CA with memory considers an extension to the standard framework of CA 
by implementing memory capabilities in cells from its own history. 

Thus to implement memory we incorporate a memory function 4>, as follow: 

^-V..,^- 1 ,^)-* (3) 
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such that t < t determines the degree of memory backwards and each cell Si G S 
being a state function of the series of states of the cell X{ with memory up to 
time-step. Finally to execute the evolution we apply the original rule as: 

V(- ■ - . «'» ■-■)-» 

Thus in CA with memory, while the mappings ip remain unaltered, historic 
memory of all past iterations is retained by featuring each cell by a summary of 
its past states from <f>. Therefore cells canalize memory to the map ip. 

As an example, we can see a case of memory <f>. The majority memory is 
defined as follow: 

4*maj * &i (4) 

where in case of a tie given by Si = So from <j> then we will take the last value 
Xi. So (f> ma j function represents the classic majority function [27] on the cells 
(xl~ T , . . . , x\) and define a temporal ring before to get finally the next 
global configuration c. 



Ahistoric t—3 t=4 t = 5 t — Q t — 8 




Figure 1 : The effect of majority memory of increasing depth on Rule 30 starting 
from a single site live cell. 

Majority memory exerts a general incrtial effect [6]. This effect, when start- 
ing from a single site live cell notably restrains the dynamics, as it is illustrated 
on Rule 30 in fig. 1. This figure shows the spatio-temporal patterns of both the 
current x state values and that of the underlying s ones. 

3 Elementary CA Rule 30 

Rule 30 was initially studied by Wolfram in [41] because of its chaotic global 
behavior; looking for a number random generator. Rule 30 is an elementary CA 
evolving in one dimension of order (2, 1). An interesting property is that Rule 
30 has a surjective relation and thus do not have Garden of Eden configurations 
[1]. In this way, any configuration has always at least one predecessor. 
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The local rule tp corresponding to Rule 30 is following: 

f 1 if 100,011,010,001 
- | o if 111,110,101,000 ' 

Generally speaking, Rule 30 displays a typical chaotic global behavior - 
Class III in Wolfram's classification. An interesting study on Rule 30 showing a 
local nested structure that repeat periodically looking invertible properties [34] . 

So, initially (fmo has a 50% of probability of states zero or one, and conse- 
quently each state appears with the same frequency. 




Figure 2: Typical behavior of Rule 30 where a single cell in state 1 leads to 
chaotic state (left). Second configuration (right) shows automaton behavior 
from a random initial condition with initial density of 50% for each state. Both 
automata evolve on a ring of 497 cells (with periodic boundary property) to 417 
generations. White cells represent state and dark cells the state 1. 

Also the evolution rule presents the following feature: if an initial configu- 
ration is covered only of state one then the configuration always evolves into 
one but if this is empty or filled with state one then this always evolves to state 
zero. Fig. 2 shows such two cases of typical evolutions in Rule 30. 

3.1 De Bruijn and subset diagrams in Rule 30 

Given a finite sequence w G S m ; such that w — w\, . . . ,w m ; let a(w) — w\, 
P(w) — u>2, • ■ ■ , w m , and ip( w ) = w\, . . . , w m _i. With these elements, we can 
specify a labeled directed graph known as de Bruijn diagram B — {N; E} asso- 
ciated with the evolution rule of the CA. The nodes of B are defined by N — S 2r 
and the set of directed edges E C S 2r x S 2r is defined as follows: 

E ={{v,w) | v,w e N, 0(v) = i>(w)}. (5) 

For every directed edge (v, w) G E, let rj(v, w) — aw G £ 2r+1 where a — a(v); 
that is, r](v, w) is a neighborhood of the automaton. In this way, the edge (v, w) 
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is labeled by ip o rj(v, w); hence every labeled path in B represents the evolution 
of the corresponding sequence specified by its nodes. Since each w £ N can be 
described by a number base k of length 2r, every node in B can be enumerated 
by a unique clement in Z fe 2r, which is useful to simplify the diagram. The de 
Bruijn diagram associated to Rule 30 is depicted in fig. 3; where black edges 
indicate the neighborhoods evolving into and the ones evolving into 1 are 
shown by gray edges. 1 



1 1 




Figure 3: De Bruijn diagram for the elementary CA Rule 30. 

Fig. 3 shows that there are four neighborhoods evolving into and four into 
1; so each state has the same probability to appear during the evolution of the 
automaton; indicating the possibility of the automaton is surjective, i.e. there 
are no Garden of Eden configurations. Classical analysis in graph theory has 
been applied over de Bruijn diagrams for studying topics such as reversibility 
[36]: cycles in the diagram indicate periodic elements in the evolution of the 
automaton if the label of the cycle corresponds to the sequence defined by its 
nodes, in periodic boundary conditions. The cycles in the de Bruijn diagram 
from fig. 3 are presented in fig. 4. 



shift in 1 gen 



Q 



Figure 4: Cycles in the de Bruijn diagram and the corresponding periodic evo- 
lution for cycle (1,2). 



1 Dc Bruijn and subset diagrams were calculated using NXLCAU21 designed by Mcintosh; 
available from http://delta.cs.cinvestav.mx/~mcintosh 
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The largest cycle in fig. 4 indicates that the undefined repetition of sequence 
Wb = 10 establishes a periodic structure without displacement in one generation 
during the evolution of Rule 30 we will than Wb is the filter 2 in Rule 30; in this 
case, the filter reported in the original Rule 30 and with memory too. Thus we 
can see how a de Bruijn diagram can recognize any periodic structures in CA 
[29, 24]. 

De Bruijn diagram is nondeterministic in the sense that a given node may 
have several output edges with the same label; a classical approach to analyse 
the diagram would be to construct the subset (or power) diagram in order to 
obtain a deterministic version for the de Bruijn diagram in the evolution rule 
[26, 38]. 




Figure 5: Subset diagram for Rule 30. 

The subset diagram is defined as S = {V, Q} where V = {P \ P C E U 0} is 
the set of nodes of S and the directed edges are defined by Q C V x V where 
for Pi, P2 E V there is a directed edge (Pi, P2) labeled by a £ S in S iff P2 is 
the maximum subset such that for every c € P2 there exists b € Pi such that 
(p(b,c) = a. 

The inclusion of the empty set assures that every edge has a well-defined 
ending node. For a CA with k states, it is fulfilled that |P| = 2 k which 
implies an exponential growth in the number of nodes in S when more states 

2 A filter is a periodic sequence that exist alone or in blocks into of the evolution, thus a 
suppression of such string produce a new view. 
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are considered. Every P £ V can be identified by a binary number showing the 
states belonging to this subset; that is, taking the states as an ordered list, the 
ones in P can be signed by a 1 and the others by 0; conforming a unique binary 
sequence identifying the subset. The decimal value of this binary number can 
be taken to get a shorter representation; where the empty set has a decimal 
number and the full subset p = X has the number 2 fe — 1. The a subset 
diagram corresponding to Rule 30 is shown in fig. 5. 

In fig. 5, the subset diagram has no path starting from the full subset (node 
15) going to the empty one (node 0). This means that every sequence can 
be produced by the evolution of the automaton and there are no Garden-of- 
Eden sequences. Thus the automaton is surjective. The subset diagram can be 
also used as a deterministic automaton for calculating ancestors of any desired 
sequences [35] by recognizing the regular expressions which may be generated 
by the corresponding automaton. Some of these expressions would be able to 
represent interesting structures as gliders [30]; however more effort is needed in 
order to get a straightforward detection of such constructions in the diagram. 

Finally, such diagrams helps to get periodic strings than eventually represent 
a general filter Wb working at the original Rule 30 and Rule 30 with memory 
(next section). Also, we will take advantage of these results to find gliders from 
these strings. 

4 Majority memory helps to discover complex 
dynamics in Rule 30 

This section reports how the majority memory <J) helps to discovers complex dy- 
namics in elementary CA by experimentation, see an introduction to elementary 
CA with memory in [8, 9, 5]. 

In this way, fig. 6 displays different scenarios where the majority memory 
works on Rule 30 and therefore extracting complex dynamics on (f) ma j- We 
should read the evolutions from left to right and up to down. Also in this case 
all evolutions use the same random initial density and the filter Wb was utilized 
in all evolutions (including the original Rule 30 evolution). 

Thus the first evolution is the original function Rule 30, i.e, without majority 
memory in fig 6. From its original evolution we can see gaps that the filter can 
clean. Traditionally was difficult distinguish such filter in Rule 30 but when 
<Pmaj w as applied then its presence was more evident. A general technique to 
get filters was developed by Wuensche in [44]. 

Initially, even values of r seems extract gliders more quickly and odd values 
fight to reach an order. Eventually the majority memory will converge to one 
stability in <j> while r increases. 

The first snapshot calculating <p ma j with t = 3 was the second evolution in 
fig 6. Nevertheless, here is not very clear how the memory could induce another 
behavior because its global behavior is similar at the original but with small 
changes. 
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Figure 6: Complex dynamics emerging in Rule 30 with majority memory <j} ma j 
from a range of values from r = 3tor = 21 although the first evolution shows 
the original function. Evolutions were calculated on a ring of 104 cells in 104 
generations with a random initial density of 50% where the initial condition is 
the same in all cases. Also the filter Wb was applied to see clearly the structures. 
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On the other hand, the third evolution with r = 4 extracts periodic pat- 
terns emerging in Rule 30 with memory. The evolution maybe does not display 
impressive gliders but it already allows to pick up more diversity in mobile lo- 
calizations on lattices of 100 x 100 cells. Thus we have enumerated and ordered 
values of r from fig. 6 based on their space-time dynamics they are responsible 
for. 

Chaotic global behavior t = 0, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21. 
Periodic patterns t = 4,6, 8, 10, 12, 14, 16, 18, 19, 20, 21. 
Collision patterns t = 6, 8, 10, 12. 

4.1 Morphological complexity in Rule 30 with memory 

Looking a way to explore global complex dynamics in Rule 30 and itself with 
memory 4> ma j- In this section we will explore some techniques. 

We evaluate morphological complexity of CA studied using morphological 
richness approach [4]. We calculate statistical morphological richness \x as fol- 
lows. Given space-time configuration of a one-dimensional CA, for each site of 
the configuration we extract its 3 x 3 cell space-time neighborhood state and 
build distribution of the space-time neighborhood states over extended period 
of the automaton's development time. 

Examples of morphological richness p are shown in fig. 7. A control case, 
when next state of a cell is calculated at random the distribution (fig. 7a) of 
space-time neighborhood states is uniform. Two-dimensional random configu- 
rations are morphologically rich. Morphology of mcmorylcss, classical, Rule 30 
automaton is characterized by few peaks in local domains distributions, where 
several space-time templates dominate in the global space-time configuration 
(fig. 7b). The statistical morphological richness p decreases. 

Incorporation of a memory in the cell-state transition rules leads to erosion 
of the distribution (fig. 7c) and thus slight increases in p. With increase of 
memory depth shape of the morphological distribution changes just slightly, up 
to minor variations in heights of major peaks (fig. 7d-f). 

Number p of 3 x 3 blocks (of states and 1) that never appear in space- 
time configurations of a cellular automaton can be used as a express estimate 
of the nominal morphological richness, the less p the more rich is a nominally 
configuration. 

The difference between statistical p and nominal p measures of morphological 
richness is that p allows to pick up most common configurations of local domains, 
while p just shows how many block of 3 x 3 states appeared in the automaton 
evolution at least once. 

For the case of random update of cell states, all blocks are present in the 
space-time configuration, p = 0. 

Memoryless automata governed by Rule 30 have p = 434 (total number of 
possible blocks is 512). When memory is incorporated in the cell-state transition 
function, first richness decreases, e.g. for r = 1, we have p = 448. Then 
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(c) 



(d) 






(e) (f) 

Figure 7: Morphological richness. CA length 1500 cells, running time 5000 
steps, (a) random update of cell states, (b) Rule 30 without memory, (c)-(f) 
with memory, (c) r = 3, (d) r = 5, (c) r = 10, (f) r = 21. 
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we observe consistent increase in complexity. Thus Rule 30 with small-depth 
memory (r = 2) p = 140, which drastically decrease to p = 68 for r = 3. The 
richness is stabilized, or rather oscillate around values p = 20 — 40 with further 
increase of memory. 

In summary, we found that majority memory increases nominal complexity 
of cellular automata but decreases statistical complexity. 

4.2 Gliders in Rule 30 with memory r = 8 

Most frequently the complex dynamics related in elementary CA is related to 
gliders, glider guns, and non-trivial reactions between localizations, for example, 
Rule 110 or 54 [28, 23]. The relevance to look for such things is the study 
of universality [14, 42]. Another additional subject was its representation by 
regular expressions and other tools [29, 24]. 




Figure 8: Gliders emerging in Rule 30 with <p ma j and r = 8. This evolution 
shows how some kinds of gliders arise and still interact from random initial con- 
ditions. The evolution was calculated on a ring of 590 cells to 320 generations, 
and an initial density of 50%. 

Among the sets of complex dynamics in Rule 30 determined by r (showed 
in fig. 6), we have chosen the memory (f> ma j with r — 8. In this way, fig 8 
illustrates an ample evolution space of its global dynamics. 

Of course, these gliders maybe are not as impressive as others from well- 
known complex rules as Rule 110, Rule 54 or other one dimensional rules [12, 
28, 23, 41, 43]. However it is interesting how 4> ma j is able to open complex 
patterns from chaotic rules. 

Nevertheless, even that Rule 30 does not offer an ample range of complex 
dynamics, it is useful to describe gliders and collisions. So, we shall illustrate 
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how a chaotic CA can be discomposed as a complex system. 




Figure 9: Set of gliders G_R3o m with memory (j> ma j and r = 8. 



Thus we classify the family of gliders and enumerate some properties. Then 
fig. 9 displays the family of gliders G#3o m = {g\, g2, 93}- Although as was 
hoped, an immediate consequence was that gliders in CA with memory have 
longer periods. 
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Table 1: Properties of gliders Gn30m with memory (f> ma j and r = 8. 

Tab. 1 summarizes the basic properties of each glider. Particularity all glid- 
ers in this domain have a constant displacement of four cells to the right and 
none glider with speed zero was founded, and yet it was more complicated to 
find a glider gun in this domain. Nevertheless, some interesting reactions were 
originated from G_R 3 o m . 

Structure w b does not have a displacement and also it is neither a glider. This 
pattern is the periodic background in Rule 30 and still represent the filter in Rule 
30. It was really hard to detect the existence of a periodic background evolving 
the original rule; but when <fi ma j was applied, a periodic pattern began to emerge 
inherited from <p_R3o- Finally this filter was confirmed with its respective de 
Bruijn and cycle diagrams (see section 3.1). 

4.3 Reactions between gliders from GRs 0rn 

Let us demonstrate simple examples of collisions between gliders. Codes for all 
gliders, necessary to generate the whole set of binary collisions, are presented 
in appendixes A and B. 

Fig. 10 shows how a stream of g\ gliders is deleted from a reaction cycle: 
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.93 + .91 .92 and g 2 + e/i -*■ g 3 - 




Figure 10: Deleting streams of 51 gliders initialized with a single 53 glider. 

To obtain such cycle by glider reaction we can code an unlimited initial 
condition as . . . g 3 ■■ -gi ■■■ -gi ■■■ -gi ■■■■ , that can be reduced as . . .53 . . . g\ ( . . . .g\ ) * (where 
a dot represent a copy of Wb)- Finally the evolution produces the cycle above, 
which is produced with nine periods of 173 and eight periods of 52- Thus each 
column presents 1135, 2270 and 3405 generations respectively. 3 . However in 
this case the codification of gliders is affected because could not natural as was 
done in other cases [28, 29, 24]. 

Since the better way to preserve ipmo and 4> ma j to code gliders as "natural," 
we will consider the codification from its original initial condition (see appendix 

3 The glider reactions were produced using OSXLCAU21 system available at http:// 
uncomp . uwe . ac . uk/genaro/OSXCASy stems . html 
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B where some strings are defined to get gliders with memory from its original 
function). 




Figure 11: Some simple reactions display how to delete (left), read (centre) and 
preserve (rigth) information with gliders in Rule 30 with memory. 

Some simple but interesting reactions from Gmom are illustrated in the 
fig. 11. First reaction is annihilation of gliders 32 and g\. Second reaction is a 
transformation g 3 glider transforms gi glider to g-i- Third reaction is a soliton- 
likc collision between gliders gi and g\. The soliton reaction between gliders is 
particularly promising because it can be used to implement computation, e.g. 
as it has been done in carry-ripple adder embedded by phase coding solitons in 
parity CA [32, 19]. 

4.4 Quasi-chemistry of gliders 

Assuming gliders g\, g 2 and g 3 are chemical species a, b and c in a well-stirred 
chemical reactor we can derive the following set of quasi-chemical reactions from 
the interactions between the gliders above: 
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a + b — > e 

a + 6 — > 2a 

, j, °- 2 
a + 6 — > c 

, °- 44 i. 
a + c — ► o 

a + cH 2 6+c (6) 

, °- 22 o 
a + c — > 3a 

0.22 

a + c — > a 

i . 0.5 , 

b+ c—> b 

b + c H 2a + b 

where reaction rates are evaluated from the frequencies of the interactions. 

We evaluated global dynamics of the quasi-chemical system (6) with constant 
volume (which reflects finite size of the automaton lattice), constant tempera- 
ture, and variable pressure using Chemical Kinetics Simulator. 41 Fig. 12 shows 
temporal dynamics of species concentrations in the system with 10 7 molecules. 

In the initial conditions when all three species are present in equal concentra- 
tions (fig. 12a) we observe exponential decay of species b and c and stabilization 
of concentration of species a. When only species b and c are present in a well- 
stirred reaction initially, they produce species c in their reactions. This leads to 
outburst in species a concentration (fig. 12b) on the background of exponential 
decline of species c and 6, until species b and c disappear and concentration of 
species a becomes constant. 

4.5 Glider machines 

Experimentally found interactions between gliders can be simplified as follows, 
depending on distance a between interacting gliders: 



(7 = 3 (7 = 4 (7 = 5 (7 = 6 

b + a^{a,b} b + a -> {0} b + a -> {0} b+a^a 
c+a—tb c + a^b c + a^>{b, c} c+a^>{b,c} 

c+b—tb c+b—tb c + b^>{a, b} c+b^>{a, b} 

Taking into account gliders' velocities (tab.l), we can construct the following 
finite state indeterministic machine with an internal state h and input state 
p, h,p € {a, 6, c, 0}. The machine can be characterised by an input-output 
transition matrix M = {mij)i,j£{ a ,b,c,%}i where for j = h f , i = p*, rriij = h t+1 . 
The matrix has the following form: 



http : //www . almaden . ibm . com/st/computational_science/ck/?cks 
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Figure 12: Dynamics of concentrations of species a (circle), b (diamond) and 
c (square) governed by reaction eqs. 6. (a) Initial concentrations of all species 
are 0.001 mole/1, (b) Initial concentrations of species b and c are 0.001 mole/1, 
species a is nil. Horizontal axes is a time in 10~ 3 sec. 
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Starting at a randomly chosen initial state and subjected to random uni- 
formly distributed input strings the machine will end in the state c with prob- 
ability \ and in the state set with probability |4. The machine starting in the 
initial h generates the string 1(h) as follows: ^(0) = 0*, 1(a) = (ab* a*)*0* ,l(b) = 
(b*a*)*®*, 1(c) = c*. 



5 Discussion 

We enriched elementary CA rule 30 with majority memory and demonstrated 
that by applying certain filtering procedures we can extract rich dynamics of 
travelling localizations, gliders, and infer sophisticated system of quasi-chemical 
reactions between the gliders. We shown that the majority memory increases 
nominal complexity but decreases statistical complexity of patterns generated 
by the CA. By applying methods of de Brujin diagrams and graph theory we 
proved that surjectivity of Rule 30 CA with memory and provide blue prints 
for future detailed analysis of glider dynamics. 




Figure 13: New family of elementary CA composed. 
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Recalling previous results on classification of one-dimensional CA [40, 43, 44] 
we envisage that introduction of majority memory 4> ma j into elementary CA 
will open a new field of research in selection of non-trivial rules of cell-state 
transitions and precise mechanics of relationships between chaotic and complex 
systems. 

This is because Rule 30 was grouped into a cluster of rules with similar be- 
haviour, that can be transformed one to another using combinations of reflexion, 
negation and complement (done by Wuensche in [40]). This way fig. 13 shows a 
diagram that explains how the original cluster for Rule 30 is presented in [40]. 
Obviously the cluster can be arbitrarily enriched using not only 4> ma j but any 
types of memory and r. Thus the dynamical complexity of automata with <j) ma j 
is the same as for the set of function {</?_R3o, (fR86, '■PRisb, ¥Ri49}, particularly 
because the local functions <^rs,q and fRU9 are responsible for leftward motion 
of gliders. 

Therefore, memory in elementary CA and others families of CA offer a new 
approach to discover complex dynamics based on gliders and non-trivial interac- 
tions between gliders. This can be substantiated by a number of different tech- 
niques, e.g. number-conservation [11, 18], exhaustive search [17], tiling [29, 22], 
de Bruijn diagrams [24], Z-parameter [44], genetic algorithms [16], mean field 
theory [25] or from a differential equations point view [13]. 
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Appendix A. Binary reactions and beyond 

We represent binary collisions in Rule 30, 4> ma j and r = 8 in the form gj space ) 
gi, where j > i, g £ Grzq and space is the interval between gliders given by the 
number of strings uib- Collisions producing e mean annihilation of gliders. Reac- 
tions are developed using increasing the distance between gliders before collision. 



Collisions of type g 2 — ► g\ 

1- 92 ^ 9i = 9i + 92 (soliton) 



o 4 

2- 92 -> .gi = e 

3- 92 -> 9i = e 



4. 92 -> 9i = 2 9i 

5- 92 9i = e 

6- 92 9i = .93 



Collisions of type gs — > g\ 
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1- .93 -» .9i = .92 

2- ff3 ^ .9i = .92 

3- 93 .9i = .92 + .93 

4- g 3 -> 5i = 5r 

5- 53 ^ 5i = 5i + 5i 

6- 53 -> 5i = .92 

7- 53 -» 5i = 5i 

8- 53 — > 5i = .92 

9- 53 — * 5i = 5i 
Collisions of type g 3 — > g 2 

1- 53 -» 52 = 52 

2- 53 .92 = 52 

3- 53 .92 = 251 + .92 

4- 53 .92 = 251 + .92 

Some other reactions with packages of gliders 

1- .92 -» 2.91 = gi 

2- 52 2ffi = 3.9i 

3. g 2 2.91 = 32 (sequence 52, 53, 52) 



4. 


6 
252 -> 5i 


= 35i 


5. 


252 5i 


= .92 


6. 


252 -> 5i 


= 51 + .93 


7. 


53 -> 2gi 


= e 


8. 


53 -> 2gi 


= 53 


9. 


53 -» 2gi 


= 51 + .92 


10. 


53 -> 2gi 


= 5? + 5i 



11- .93 -> 2gi = g\ 

12. g 3 ^ 2 9l = 2 9l 

13. 2g 3 — > g 2 = 32 (wall) 

Appendix B. Coding gliders G^om 

We can enumerate strings conforming gliders in Rule 30 with 4> ma j and t = 8, in 
given initial conditions and using "phases" (we omit strings that do not produce 
gliders). 

Note that such strings evolve initially with (PR30 and a value of r given, 
then 4> ma j will open these strings when memory works. Thus we can code 
initial conditions with gliders in CA with memory, that also was implemented 
in OSXLCAU21 system to get our simulations. 



5i glider 


52 glider 


33 glider 


1-51 

111-51 
10000 - 5i 
11110 - 5i 
100 - 52 
11110 - 5i 
10000 - 5i 
11001 -52 + 53 


100 - 52 

111-51 
10000 - 5i 
11110 - 5i 

1100110 - 52 +53 
1010110 - 5i 
1101100 - 52 


110110 - 251 join 
101110 - 53 
111100 - 2 5 i 
100011 - 2 5 i join 
IOII-53 



Table 2: Strings evolving in gliders of Git30m- 



Tab. 2 enumerates strings for each glider represented as a tiling so we know 
their "phases" [28, 29]. Although in this case, it was difficult to classify such 
strings as regular expressions because not all strings from the tiling representa- 
tion evolve into gliders. 
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